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ABSTRACT 

We use the Arepo moving mesh code to simulate the evolution of molecular clouds 
exposed to a harsh environment similar to that found in the galactic center (GC), in 
an effort to understand why the star formation efficiency (SFE) of clouds in this envi¬ 
ronment is so small. Our simulations include a simplified treatment of time-dependent 
chemistry and account for the highly non-isothermal nature of the gas and the dust. 
We model clouds with a total mass of 1.3 x 10^ Mq and explore the effects of varying 
the mean cloud density and the virial parameter, a = F^kin/|^pot|- We vary the latter 
from G = 0.5to(a = 8.0, and so many of the clouds that we simulate are gravitation¬ 
ally unbound. We expose our model clouds to an interstellar radiation field (ISRF) 
and cosmic ray fiux (CRF) that are both a factor of 1000 higher than the values found 
in the solar neighbourhood. As a reference, we also run simulations with local solar 
neighbourhood values of the ISRF and the CRF in order to better constrain the ef¬ 
fects of the extreme conditions in the GC on the SFE. Despite the harsh environment 
and the large turbulent velocity dispersions adopted, we find that all of the simulated 
clouds form stars within less than a gravitational free-fall time. Increasing the virial 
parameter from a = 0.5 to g = 8.0 decreases the SFE by a factor ^4 — 10, while 
increasing the ISRF/CRF by a factor of 1000 decreases the SFE again by a factor 
^2 — 6. However, even in our most unbound clouds, the SFE remains higher than 
that inferred for real GC clouds. We therefore conclude that high levels of turbulence 
and strong external heating are not enough by themselves to lead to a persistently low 
SFE at the center of the Galaxy. 
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1 INTRODUCTION 


Milky Way and also in distant galaxies (Kennicutt & Evans 


Understanding star formation is an important task in the¬ 
oretical astrophysics (see e.g., Mac Low Klessen||2004| 


Scalp Elmegre^|2004| [Elmegreen Scalo||2004| |McKee| 

fc Ostriker|20Q7 Ballesteros-Paredes et al.|2007 ). Stars form 
in molecular clouds (MC) in the interstellar medium (ISM) 
due to gravitational contraction of overdense gas regions (see 


also the lecture notes by Klessen & Glover 2014). Several 
studies in the past have revealed that the number of stars 
formed per unit time (referred to as the star formation rate 
or SFR) is proportional to the amount of gas in the star¬ 
forming region (see e.g., Schmidt|I959 Kennicutt 1 1998 ) and 
that the SFR surface density shows a power-law dependence 
on the gas surface density. This is known as the Schmidt- 
Kennicutt relation and is an important empirical result in 
star formation theory. It appears to hold in the disk of our 


2012). However, there is still some debate about whether a 


universal relationship holds for all galaxies or not (see, e.g. 
Shetty et al.polslp^ l. 


As shown by Longmore et al.| (2013a), the inner 500 pc 
region of our Galaxy, known as the Central Molecular Zone 
(CMZ), contains the largest reservoir of dense gas in the 
Milky Way, with densities of order of several 10^ cm“^ and 
a total mass of about 10^ M©. This would be enough to 
form several Orion-like clusters, but the measured SFR is 
significantly smaller ( Longmore et al^ 2013a). This is a di¬ 
rect contradiction to the Schmidt-Kennicutt relation, which 
would predict a high SFR in the center of our Galaxy. The 
key questions are, why on one hand the observed SFR is 
suppressed by a factor of > 10 ( Longmore et al.] 20I3a) and 
on the other hand which physical processes regulate star 
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formation in the central few hundred parsecs of the Milky 


Way (see, e.g. Longmore et ah 2013a, Kruijssen et ah 2014, 
Krumholz & Kruijssen 2015). 

Several studies have tried to solve this problem of star 
formation in such an extreme environment, which challenges 
current star formation theories. For example, |Kruijssen et al.| 
(2014) studied the impact of several mechanisms on the SFR 


on different physical scales, e.g. the very strong radiation 
field, magnetic fields, turbulent pressure, galactic tides or 
feedback. They argue that star formation could be episodic 
due to a gradual build-up of dense gas by spiral instabilities 
or that variations in the rates of gas flows into the CMZ 
might significantly alter the star formation process. |Long-| 
more et al.| (2013a) state, for example, that the order of 
magnitude higher internal velocity dispersion could disrupt 
compact regions before they are able to go into gravitational 
collapse. However, the question of what physical processes 
are most important for regulating the SFR in the galactic 
center (GC) remains unresolved. 

One important way in which we can distinguish between 
different models for the regulation of star formation in the 
GC is by examining whether they produce low star forma¬ 
tion efficiencies by reducing the formation of stars in the in¬ 
dividual dense clouds present in the CMZ, or whether they 
instead merely regulate the rate at which molecular gas is 
allowed to form these dense clouds. The presence of a num¬ 
ber of massive, dense clouds in the CMZ that do not appear 
to be currently forming stars provides a hint that the star 
formation efficiencies of individual clouds in the CMZ might 
be low. An example of such a cloud, which is almost devoid 
of star formation, is also known as “The Brick” (see e.g. 


Glisten et al.||1981||Lis et al.||1994||Lis & Menten||1998||Lis 

et al.pOOl 

Molinari et al.||2011| jlmmer et al.||2012 

Long- 

more et al. 

2012||Kauffmann et al.||2013| Clark et a 

. 2013[ 

Johnston et al.|2014| Rathborne et al.|2014||Kruijssen et al.j 

2015||Pillai 

et al.|2015| as well as [Longmore et al. 2013b for 


a discussion of several other dense starless clouds close to 
the Brick that have similar properties). On the other hand, 
it is also possible that we are just observing the Brick and its 
neighbouring starless clouds at a very early time in their evo¬ 
lution, before they have started to form stars. This question 
is difficult to resolve observationally, but numerical simula¬ 
tions can help us to understand which of the possibilities is 
more likely. 

In this paper we investigate the impact of a strong in¬ 
terstellar radiation field (ISRF) and a high cosmic ray flux 
(CRF) on the SFR in clouds with different initial number 
densities and varying levels of turbulence. We change the 
amount of turbulent kinetic energy F^kin with respect to the 
potential energy F^pot, as parameterized by a = F^kin/|F^pot | 
(see Section |2.3| ) . We focus on the question of whether such 
a high ISRF/CRF combined with a high level of turbulence 
could be the main physical drivers to suppress star forma¬ 
tion in the GC. Therefore, we adopt environmental condi¬ 
tions that are similar to those experienced by a typical GC 
cloud. In particular, we adopted values for the ISRF strength 
and the CRF comparable to those inferred for the Brick. We 
model the behaviour of different MCs exposed to this envi¬ 


ronment using the Arepo moving mesh code (Springel 2010) 


and explore the effect of changing their density and turbu¬ 
lent velocity dispersion. For technical reasons, the clouds we 
model have lower densities than typical GC clouds, and so 


should be less likely to form stars than real GC clouds. De¬ 
spite this and despite the harshness of their environment, it 
proves to be very difficult to suppress star formation within 
the clouds, as we will see in the following sections. 

The structure of our paper is as follows. In Section 
we present the simulations and the methods used in this 
paper. In Section we present the results of our studies 
with different virial a parameters and various initial number 
densities. We discuss our results in Sectionand present our 
conclusions in Section]^ 


2 METHODS AND SIMULATIONS 


2.1 Computational method 

Our simulations are performed using the moving mesh code 
Arepo (Springel 2010), which uses an unstructured mesh 


defined by the Voronoi tessellation of a set of discrete points. 
We make use of a detailed atomic and molecular cooling 


function, described in detail in Glover et al. (2010) and 


Glover & Clark (012b), and a simplified treatment of the 


molecular chemistry of the gas. Our chemical treatment is 


based on the work of Nelson & Langer (1997) and Glover 


& Mac Low (2007), and allows us to follow the forma¬ 


tion and destruction of H 2 and CO self-consistently within 
our simulations. Full details of the chemical model with a 
description of how the chemistry interacts with the ISRF 
via the TreeCol algorithm can be found in [Clark et al.| 
(2012). Examples of the use of our chemical model with the 
Arepo code can be found in Smith et al.|(2014a) and Smith 
eTaL] (2014b). 


We assume that the gas has a uniform solar metallic- 
ity and adopt the standard ratio of helium to hydrogen, and 
abund ances of carbon and oxygen taken from jSembach et aiL\ 
(2000), i.e. xc = 1.4 x 10“^ and xo = 3.2 x 10“ , where xc 


and Xo are the fractional abundances by number of carbon 
and oxygen relative to hydrogen. However, we have to keep 
in mind that the CMZ has actually a super-solar metallicity. 
Nevertheless, we use a uniform solar value in order to be con¬ 
servative regarding the cooling and star formation rates in 
our runs. At the start of the simulations, hydrogen, helium 
and oxygen are in atomic form, while carbon is assumed to 
be in singly ionized form, as C^. We also adopt the stan¬ 


dard local value for the dust-to-gas ratio of 1:100 (Glover 


et al. 2010), and assume that the dust properties do not 


vary with the gas density. The cosmic ray ioniz ation rate of 
atomic hydrogen is set to C = 3 x 10“^^ s“^ (Clark et al. 
2013), which is a factor of ~ 1000 higher than the value in 


the solar neighbourhood ( Yusef-Zadeh et ar||2007 ). For the 
incident ultraviolet radiation field, we adopt the same spec¬ 


tral slope as given in Draine (1978). We denote the strength 


of the Draine ISRF as Co = 1 and perform simulations with 
a field strength Go = 1000 ( Clark et al.||2013 ). The Draine 
field has a strength Go = 1.7 in Habing (1968) units, cor- 


FT- 


Furthermore, as a reference, we also run simulations with 
local solar neighbourhood values of the ISRF and the CRF 
in order to explore the effect of a different radiation field on 
the SFE. In this case, we set the field strength of the ISRE 
to Co = 1 and the cosmic ray ionization rate of atomic hy¬ 
drogen to C = 3 X 10“ 
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refinement criterion, which is active over the whole simula¬ 
tion period in order to accurately refine dense and collapsed 
gas regions in the box. We use a constant number of 8 cells 
per Jeans length, which is sufficient to avoid artificial frag¬ 
mentation ( Truelove et al.||1998| Greif et ffi^|2Qll ). 


2.2 Sink particles 


Furthermore, we make use of a sink particle implementation 
( Greif et al.||2011 ) based on the prescription in Bate et ah 
(1995) and Jappsen et ah (2005), to track the star formation 


process during the simulations. Before a cell is turned into a 
sink particle, it undergoes a series of tests. First, the particle 
must reach a critical density threshold of nthresh ~ 10^ cm“^. 
This value was chosen to be higher than the typical post¬ 
shock densities found in the clouds to ensure that we do not 
attempt to form sinks in regions that are not gravitationally 
collapsing. For our more extreme run with no = 10^ cm“^ 
we adopt a higher threshold of nthresh ~ 10^ cm“^ in order 
to account for the two orders of magnitude higher density 
on average. The second test is to check whether the new sink 
particle is sufficiently far away from any other sink particle, 
measured in terms of one accretion radius race- The third is 
to check whether the size of the cell is less than the accretion 
radius of the sink particle that it will become. 

We set the accretion radius to a constant value of 
Tacc ~ 0.01 pc, roughly corresponding to the scale of the 
thermal Jeans length for a sink particle formation threshold 
of nthresh ~ 10^ cm“^ at a mean gas temperature of ~ 20 K. 
For our model with a threshold of nthresh ~ 10®cm“^, we 
adopt the same value for the accretion radius. This is also a 
reasonable estimate for this higher density model, since more 
checks will guarantee that the gas is bound and collapsing 
onto the sink particle. In addition, we explore whether a dif¬ 
ferent value of race might have a significant impact on the 
SFEs by running a simulation in which we increased race by 
a factor of 10. We find that in this case, the SFE increases 
by around a factor of 1.1 — 1.3. In view of the various error 
sources for estimating a star formation efficiency per free-fall 


time (see also the discussion in Section 3.1), we find that a 
different accretion radius does not significantly affect the re¬ 
sults of this paper. This is because we are mainly interested 
in estimating the total mass that goes into gravitational col¬ 
lapse instead of measuring precise SEEs. 

Once these preliminary criteria are fulfilled, more checks 
guarantee that this potential sink particle is in a correct dy¬ 
namical state. In a first test, we require that the particle is 
sub-virial, i.e. it needs to fulfill the condition a ^ 0.5, where 
a is the ratio of kinetic and gravitational energy. Second, we 
ensure that a-\- ^ ^ 1, where /3 is the ratio of rotational and 
gravitational energy. Third, we require div(a) < 0, where a 
is the acceleration, which ensures that the particle is not 
tidally disrupted or bouncing. If all these conditions are 
achieved, the local gas condensation can become a sink. 


2.3 Important quantities 


In this study, we use the virial a parameter to regulate the 
amount of turbulent kinetic energy in the simulation do¬ 
main. We define this as 


_ -Ekin 

I Epot I 


( 1 ) 


where Ekin and Epot denote the total kinetic and potential 
energy in the box at the start of the simulation, given via 


and 


1 2 
Ekin — ^MioiCTy 


J^pot - 


5E 


( 2 ) 

(3) 


In this context, ay denotes the turbulent 3D velocity disper¬ 
sion, Mtot and R the total mass and radius of the (initially) 
uniform sphere and G the gravitational constant. Eurther- 
more, we use 


Tftot — 


4 

3 




(4) 


where p is the initial mass density. However, we note that the 
virial parameter is also often defined via a — 2Ekin / 1 Epot I 
in the literature, which is different from the notation used 
above. In our definition, a value of ct = 0.5 defines virialized 
clouds, ct = 1.0 denotes clouds with energy equipartition 
and ct > 1.0 describes clouds that are (highly) unbound. 
In addition, we can estimate the crossing time from these 
quantities via Gross ~ Rjerv 

Eurthermore, we quantify the amount of gas mass being 
converted to stars (i.e. sink particles) due to gravitational 
collapse within one free-fall time tff as 


Cff = tff • 


Mtot 


(5) 


where M* is the star formation rate in the computational do¬ 
main averaged over tff ( Krumholz McKee|200'5 Krumholz 


fc Tan|20Q7 Murray|2011 ). The parameter is thus a mea¬ 

sure of the star formation efficiency (SEE) of each model. 
The free-fall timescale itself is defined via 


In this paper, we compute Cff for all clouds individually and 
compare our results to the average SEE inferred for the 
galactic center region. It is therefore important to emphasize 
that our study is aimed at testing the idea that the SEE at 
the galactic center is low because the SEE of the individ¬ 
ual dense clouds is low. With our isolated cloud models, we 
cannot test the competing idea that the SEE of the region 
is low because it is difficult to form dense clouds there. 


2.4 Initial conditions and model parameters 

Eor simplicity, we assume that the cloud is initially spherical 
and embedded in a low-density environment. Eurthermore, 
we apply periodic boundary conditions. We perform two sets 
of runs in which we vary the virial a parameter, using values 
a = 0.5,1.O, 2.0,4.0 and 8.0. In the first set of runs, we take 
an initial hydrogen nuclei number density of no = 100 cm“^ 
for the cloud. In the second set of runs, we take a number 
density of 1000 cm“^. In all of our simulations, the density 
of the gas surrounding the cloud is I cm“^. We use a total 
mass of Mtot = 1-3 x 10^ M© and an initial total number of 
2 X 10® cells for the whole box, including the low-density 
regions. The initial cell mass within the uniform cloud cor¬ 
responds to ~ 2 M©. 

Given these parameters, we can compute the radii and 
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Model name 

Virial a 

Initial density no 

Radius R 

Velocity ay 

Free-fall time tfj 

Crossing time fcross 

ISRF & GRF 



[cm-3] 

[pc] 

[km/s] 

[Myr] 

[Myr] 


GG-0.5-100 

0.5 

100 

19.1 

3.7 

4.4 

5.1 

1000 

GC-l.0-100 

1.0 

100 

19.1 

5.2 

4.4 

3.4 

1000 

GG-2.0-100 

2.0 

100 

19.1 

7.4 

4.4 

2.5 

1000 

GG-4.0-100 

4.0 

100 

19.1 

10.5 

4.4 

1.8 

1000 

GG-8.0-100 

8.0 

100 

19.1 

14.7 

4.4 

1.3 

1000 

GC-O.5-1000 

0.5 

1000 

8.9 

5.4 

1.4 

1.6 

1000 

GC-l.0-1000 

1.0 

1000 

8.9 

7.6 

1.4 

1.1 

1000 

GC-2.0-1000 

2.0 

1000 

8.9 

10.8 

1.4 

0.8 

1000 

GC-4.0-1000 

4.0 

1000 

8.9 

15.3 

1.4 

0.6 

1000 

GC-8.0-1000 

8.0 

1000 

8.9 

21.6 

1.4 

0.4 

1000 

GG-16.0-10000 

16.0 

10000 

4.1 

44.9 

0.4 

0.09 

1000 

SOL-O.5-1000 

0.5 

1000 

8.9 

5.4 

1.4 

1.6 

1 

SOL-8.0-1000 

8.0 

1000 

8.9 

21.6 

1.4 

0.4 

1 

The Brick 

1.0 

7.3 X 10-^ 

2.8 

16.0 

0.34 

0.17 

- 100 - 1000 


Table 1. Initial conditions for our different cloud models. For each run, we list the virial ol parameter, the initial number density, the 
spherical cloud radius, the velocity dispersion, the free-fall time, the crossing time and the scaling factor of the interstellar radiation field 
(ISRF) and the cosmic ray fiux (CRF) relative to the solar neighourhood value. For comparison, we also list the parameters for the GC 
cloud GO.253+0.016, known as “the Brick”, which are given in Table 2 in |Longmore et n7| ( |2012| . 


the initial velocity dispersions of the clouds. The cubic side 
lengths of the total simulation domain are set to 5 x the indi¬ 
vidual cloud radii. All clouds are initially at rest and placed 
in the center of the box. We use an initial random velocity 
field with a power spectrum of P(/c) oc that consists 
of a natural mixture of solenoidal and compressive modes, 
which decays throughout the simulation. The timesteps be¬ 
tween the individual snapshots are Atgnap ~ 22kyr and 
Atsnap ~ 7kyr for models with initial number densities of 
no = 100cm“^ and 1000cm“^, respectively, corresponding 
to ~ 200 snapshots per simulation in total. 

In addition, we have also modelled an even more ex¬ 
treme cloud with an initial number density of no = 10^ cm“^ 
and a virial a parameter of a = 16.0. In this case, we set 
Atgnap ~ 2kyr for the same total number of time snapshots. 
However, running further simulations with no ^ 10^ cm“^ 
and OL < 16.0 is computationally prohibitive, owing to the 
rapid rate at which stars form in these models. For this 
reason we only focus on this one more extreme numerical 
model. Nevertheless, as we discuss in more detail later, such 
a high OL run alone is still enough for our study, since even in 
this extreme case, the cloud still forms a substantial number 
of stars. Furthermore, regarding the densities in the CMZ, 
we have to keep in mind that a significant number of GC 
clouds have densities higher than those modelled in our sim¬ 
ulations, which would lead to even larger turbulent velocities 
in our simulations for the same constant total mass. Table [T] 
gives an overview about the different numerical models and 
their initial conditions. 


cf 


100 


10 


0.1 




•^O.seedO 
*^ 0,86661 
'^0,seed0 = 


. 


: 100 cm 
: 100 cm 


1000 cm 


Eff ~ exp(-0!20 X a) 
Eff ~ exp(-0.21 X a) 
~ exp(-0.36 X a) 

X 


.X— 

-e- 


0 1 


7 8 9 10 


Virial parameter a = / |EpQt| 


Figure 2. Star formation efficiencies per free-fall time against 
the corresponding virial a parameter for our fiducial models with 
no,seedO = 100 cm“^ and no,geedO = 1000 cm“^ (Table and 
for one model using a different seed o f th e turbulent velocity 
field with u-o^geedi = 100 cm“^ (Table [At] >. We generally find 
a decreasing trend of the SFEs with higher a for all density 
models. Although there are slight differences in the individual 
SFEs between the two models with uo,geedo = 100 cm ^ and 
no,geedi = 100 cm“^, we nevertheless find the same general trends 
of decreasing Cff with increasing a. To illustrate the different 
trends, we also fit exponential functions oc exp(—ca) to the mod¬ 
els and list the corresponding slopes in the plot. Error bars in¬ 
dicate variations of the SFE of ~ 30%, which we conservatively 
estimate in Section lO 


3 RESULTS 

3.1 Analysis of the model clouds 

Fig.[T] shows logarithmic column density maps of the dif¬ 
ferent model clouds presented in Table As expected, the 


shape of all the clouds strongly depends on the specific value 
of a in each simulation. Tablegives an overview of the star 
formation efficiencies and the number of sink particles being 
formed in each model cloud with varying o, also illustrated 
in Fig. In the Figure, we also show the star formation 
efficiencies that we derive from simulations performed using 
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Figure 1. Logarithmic column density maps in z-direction for all cloud models with initial number densities of no = 100 cm“^ (left 
column) and 1000cm“^ (right column) for different virial parameters ol = 0.5,1.0, 2.0, 4.0, 8.0 (from top to bottom), showing an extract 
of the central cloud regions. In each plot we give the number of sink particles, the estimated star formation efficiency per free-fall time 
and the simulation time. Sink particles are formed during the simulations and marked with white dots in each map. The side length of 
each box shown above corresponds to 74.5 pc and 44.5 pc for the 100 cm“^ and 1000 cm“^ models, respectively. 
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Figure 3. Left column: Mass-weighted volume density PDFs at different times for our five virial a parameters, a = 0.5,1.0, 2.0, 4.0 and 
8.0. Right column: Same as left column, but with column density PDFs. In order to compute the column density PDFs, our simulations 
were projected onto a regular 1024^ map. All plots are computed for our models using an initial number density of no = 100 cm“^. 
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Model name 

^At 

eff 

A^sink 

u 

lend 

At 


[%] 

[%i 


[Myr] 

[Myr] 

[Myr] 

GC-0.5-I00 

8.2 

25.7 

3473 

1.30 

2.70 

1.40 

GC-l.0-100 

8.5 

28.5 

4699 

1.80 

3.12 

I.3I 

GC-2.0-100 

8.5 

I9.I 

4166 

1.84 

3.80 

1.96 

GG-4.0-100 

6.2 

14.8 

2511 

2.56 

4.44 

1.84 

GC-8.0-100 

2.3 

6.1 

644 

2.75 

4.44 

1.65 

GG-0.5-1000 

4.5 

10.9 

2966 

0.42 

I.OO 

0.58 

GC-1.0-1000 

4.6 

8.0 

2531 

0.28 

1.10 

0.82 

GC-2.0-1000 

6.8 

8.0 

3218 

0.22 

1.40 

1.18 

GC-4.0-1000 

2.6 

2.9 

1035 

0.16 

1.40 

1.24 

GC-8.0-1000 

0.6 

0.7 

180 

0.12 

1.40 

1.28 

GC-I6.0-I0000 

0.7 

0.7 

1406 

0.02 

0.40 

0.38 

SOL-O.5-1000 

5.1 

13.2 

16248 

0.21 

0.75 

0.54 

SOL-8.0-1000 

4.3 

4.6 

8617 

0.08 

1.40 

1.32 


Table 2. The star formation efficiencies e^t give the amount of 
gas being converted to sink particles within a time At = tend “ 
t* < tff, where tend denotes the end of our simulation and t* the 
time when the first star forms. A^gink gives the total number of 
sink particles formed during the time At. To calculate the SFEs 
per free-fall time, Cff, we extrapolate based on CAt? assuming that 
the SFR between tff and tend is the same as between and tend- 
Some simulations with lower a values are not evolved until one 
free-fall time, which is due to the high computational costs of the 
individual runs. Nevertheless, we let all those simulations evolve 
until a CAt of at least ~ 4% is reached. 


the same values of no and a, but with a different random 
seed for the turbulent velocity field, as explained in more 
detail in Appendix]^ Some of the simulations with less tur¬ 
bulent kinetic energy (i.e. those with lower a values) were 
not evolved until the end of one free-fall time. This is be¬ 
cause gravitational collapse and star formation become very 
efficient in these simulations, driving up the computational 
cost due to the need to refine many high density regions. 
In these rapidly star-forming runs, we stop our simulations 
when they either reach one free-fall time or when the compu¬ 
tational cost of continuing becomes excessive. In most cases, 
this occurs once ~ 8% of the gas has formed stars, but in 
runs GC-0.5-1000 and GC-1.0-1000 the computational cost 
climbs so steeply as the gas collapses that we are forced to 
stop when only ~ 4% of the gas has formed stars. We note, 
however, that we do not include the effects of feedback from 
young stars (see Section 4.4), and so the star formation rate 
is likely to be overestimated at late times in all of our runs. 

We also have to keep in mind that we start with an ide¬ 
alized, spherical and uniform gas distribution at the begin¬ 
ning, rather than with an MG that is already in an evolved 
physical state. Measuring the star formation efficiencies from 
the beginning of the simulation is therefore a questionable 
procedure, since our results would be strongly affected by 
the initial geometry of the MG. Thus, we give two different 
star formation efficiencies in Table The first of these, CAt, 
denotes the fraction of gas that is converted to stars within 
a time interval At = tend — t*, where tend is the time when 
we stop the simulation and t* the time when the first star 
forms. The other, Cff, is the SFE per free-fall time, which is 
computed by extrapolating CAt to one free-fall time. That 


means we evaluate 


eff = CAt 


tff 


(7) 


We also note that the number of sink particles A^sink that 
form in each simulation might depend on the specific choice 
of the sink particle formation threshold nthresh • This number 
therefore has to be treated with caution, since it is proba¬ 
bly not converged, given our resolution, while the SFR is 


converged (see also Glover & Glark 2012a). However, we are 
not aiming at resolving the IMF in this study, but instead 
want to obtain the total mass that goes into gravitational 
collapse, which is correctly described given our numerical 
setup. 

Furthermore, we note that the largest source of error in 
estimating eff is caused by the extrapolation of CAt to eff. In 
general, we assume that the mass accretion rates stay con¬ 
stant through the remaining time of extrapolation. Strictly 
speaking, this is true only for our a = 8.0 and a = 4.0 
models, as we will see later in Section [3.3| Further smaller 
errors in estimating converged values of eff are caused by 
the concrete realization of the turbulent velocity field (see 
also Appendix 0- as well as by the specific choice of the 
sink particle threshold and the accretion radius. Altogether, 
we think that an error of ~ 20 — 30% is a conservative es¬ 
timate. This is acceptable, because we are primarily inter¬ 
ested in analyzing how different levels of turbulence, density 
and the ISRF/GRF affect the formation of stars instead of 
measuring exact and converged values of the individual star 
formation efficiencies. 

In general, we find active star formation for all three ini¬ 
tial densities and for all virial a parameters, even in clouds 
which are unbound due to a high value of a. The highest 
efficiencies are found in the models with small virial param¬ 
eters, which have the lowest turbulent velocity dispersion of 
all our models (see Table [^. The lowest efficiencies are ob¬ 
tained in the GG models with large virial parameters, which 
are runs with the highest turbulent velocity dispersions. Re¬ 
garding our extreme GG-I6.0-I0000 model, we find a SFE 
of eff ~ 0.7% even in this case, although the internal veloc¬ 
ity dispersion is very high. Moreover, we generally observe 
a decreasing star formation rate per free-fall time eff with 
increasing a (see Fig.[^, which we would expect as well due 
to the increasing amount of turbulent kinetic energy in the 
simulation domain. 

We quantify this and fit an exponential law to the data 
points shown in Fig. [fusing a x^-fit. Our fitting function is 
defined via 


eff oc exp(—ca), (8) 

where a is the virial parameter and c is a constant, depend¬ 
ing on the model. For our model with no = 100 cm“^, we 
find c ~ 0.20 ± 0.02, while in the n = 1000 cm“^ model we 
obtain a steeper slope c ~ 0.36 ± 0.07 (see Fig.[^. Thus, we 
find strong evidence that eff depends not only on the virial 
state of the cloud, but also on its density. Gomparing the 
values for the SFE for our two density models, we find that 
eff ~ which means that the SFE per free-fall time 

is smaller in high density clouds than in clouds with lower 
density. Interestingly, this is the same scaling with density 
as the free-fall time itself, which suggests that the change in 
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Cff is driven largely by the change in the free-fall time, rather 
than by a systematic change in the star formation rate. 

In addition, we also run two simulations, SOL-0.5-1000 
and SOL-8.0-1000, adopting the lower solar neighbourhood 
values for the ISRF and the CRF, in order to compare the 
SFEs we measure in this quiescent environment to those 
measured for the much harsher GC environment. For the 
virialized cloud (SOL-0.5-1000), we hnd an SFE per free-fall 
time that is around 20% larger than in our corresponding 
GC model (GC-0.5-1000), demonstrating that for gravita¬ 
tionally bound clouds, the much stronger heating present in 
the GC has little effect on the star formation rate. It is no¬ 
table, however, that we form far fewer sink particles in our 
GC run than in the corresponding solar neighbourhood run, 
suggesting that the sinks that do form must be systemati¬ 
cally larger. Whether this also leads to a systematic change 
in the initial mass function of the stars forming in this envi¬ 
ronment remains to be seen; unfortunately, our resolution is 
too low to allow us to properly address this question. In our 
runs with a — 8.0, we hnd a much larger difference between 
the solar neighbourhood and GC runs. In the solar neigh¬ 
bourhood run, increasing a from 0.5 to 8.0 decreases by 
less than a factor of three, whereas in the corresponding GC 
runs, the change in Cff is closer to a factor of ~ 16. There¬ 
fore, the combination of high turbulent velocities and strong 
heating is much more effective at suppressing star formation 
than either effect individually. 


3.2 Analysis of the volume, temperature and 
column density Probability Distribution 
Function (PDF) 

Fig. and show the mass-weighted volume density and 
the column density PDF for both density models with 
no = 100 cm“^ and no = 1000 cm“^ at different times in 
the cloud evolution for all virial a parameters. In order to 
compute the column density PDF, our simulations were pro¬ 
jected onto a regular 1024^ map. We have chosen various 
time snapshots in the evolution, so that all PDFs properly 
reflect the different physical states of the cloud. At later 
times, however, we have removed some of the low a runs 
from the plots, since these runs are not evolved until the 
one free-fall time limit. Nevertheless, we can still compare 
the more evolved PDFs to the PDFs of the low a runs from 
the snapshot shown above at an earlier time. All plots only 
include the remaining total gas without the mass already 
converted to sink particles. 

The mass-weighted volume density PDF shows two pro¬ 
nounced peaks in each of the plots, suggestive of a two-phase 
medium. The peaks are found at densities oi n ^ 1 cm“^ 
and n ^ 10^ cm“^ in the no = 100 cm“^ model and at 
n ~ 10cm“^ and n ~ 10^ cm“^ in the no = 1000 cm“^ 
model, respectively. These peaks in the bimodal PDF can 
be referred to as those of the diffuse ISM and of the denser 
regions of the MG, which make up most of the mass. We 
also note that both peaks are well below the density thresh¬ 
old for sink particle formation. Moreover, both the density 
variance and the mean density generally increase as the sim¬ 
ulations evolve with time in our models. In addition. Fig. 
shows the corresponding mass-weighted temperature PDF 
for both density models. We find that most of the gas has 
temperatures of several 100 K for all virial parameters, ow- 




Figure 6. Mass of gas M* converted to stars (sink particles) as 
a function of time for the different virial a parameters for our 
models with initial number densities of no = 100 cm“^ (top) and 
no = 1000cm“^ (bottom). 


ing to the strong heating by the ISRF and the GRF. This 
stands in contrast to the situation in local MGs, where a 
large fraction of the gas mass has T ~ 10 — 20K. 

Analysing the different temporal evolutions of the PDF 
for the various a values, we find that the high density regions 
of runs with larger a values are denser than runs with lower 
a values during the first ~ 25% of the free-fall time. This 
is due to the higher internal velocity dispersion, which can 
compress the gas more effectively up to higher densities at 
the beginning of our runs. At later times, the high density 
regions of simulations with lower a values become denser 
due to gravitational contraction of the medium. Gomparing 
the PDF of the two density models in Fig.[^and|^in general, 
we hnd that both show a similar shape, except a signihcant 
shift to larger density values for the 1000 cm“^ simulations 
relative to the 100 cm“^ simulations, which is due to the one 
order of magnitude difference in the initial number density. 
We also note that we observe a power-law tail at the higher 
end of the different column density PDF, once star formation 
has set in, consistent with previous work in this field, e.g. 


Klessenj (2000a), Kritsuk et al. (2011), Federrath & Klessen 
( 2013| ) , [Schneider et al. (2013) or Rathborne et al. (2014). 
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Figure 7. Mass accretion rates dM/dt for the corresponding den¬ 
sity models and different virial a parameters given in Fig. 


3.3 Sink particle formation 

Fig. 1^ shows the mass of gas that is converted to stars (sink 
particles) as a function of time for the different a parameters 
and densities. Depending on the amount of kinetic turbulent 
energy in the box, we find differences in the temporal evo¬ 
lution of star formation. In our models with a lower initial 
number density of no = 100 cm“^, we find star formation 
to be triggered by the global collapse of the MC due to 
gravitational compression. This is more effective for clouds 
with lower virial parameter. Hence, in this case, turbulence 
delays and suppresses star formation. For our other mod¬ 
els with higher initial number densities of no = 1000 cm“^, 
we find star formation to be triggered by local compression 
of the gas due to highly turbulent motions, leading to the 
rapid formation of stars in models with high a values. Over¬ 
all, star formation is suppressed at high <a, but not delayed 
as observed in our low-density model. Instead we note that 
turbulence can actually trigger and accelerate star formation 
in localized patches of the cloud (see also the discussion by 
Mac Low Kless^|20Q4 ). 

Sink particles can form once all formation criteria pre¬ 
sented in Section [2^ are fulfilled. However, highly turbulent 
motions can locally compress the gas above our sink parti¬ 
cle formation threshold more quickly at higher values of the 
density and the virial parameter. Hence, under our assump¬ 
tions of sink particle formation, we find that star formation 
can be triggered by different physical processes, i.e. either 


by shock compression of the gas leading to overdense regions 
which fulfill all sink criteria, or by global gravitational col¬ 
lapse. In the latter case, the global collapse dominates the 
internal velocity dispersions of the cloud, which is e.g. the 
case in quiescent MCs in the Galaxy. However, we assume 
that in a CMZ-like environment with highly turbulent mo¬ 
tions, star formation is mainly triggered by turbulent shock 
compression rather than by global gravitational collapse of 
the cloud as a whole. 

Fig.j^shows mass accretion rates for the corresponding 
density models and different virial parameters of our simula¬ 
tions. Depending on the amount of turbulent kinetic energy, 
we find a strong dependence of the accretion rates on the 
a parameter. As expected, the gas can be accreted more 
effectively in case of low a values, i.e. when the kinetic en¬ 
ergy is low and the cloud tends to be more gravitationally 
bound, which holds for all our density models. Furthermore, 
for those models with a low value of the virial parameter, 
we find that star formation continues to accelerate during 
the whole run. If a increases, star formation gets less ef¬ 
ficient, the mass accretion rates drop as well and become 
approximately constant at later times. 


4 DISCUSSION 


4.1 Comparison to star formation in the GC 

In this study, we have adopted environmental conditions 
similar to those experienced by a typical GC cloud in order 
to see whether high levels of turbulence as well as a much 
stronger ISRF and CRF will lead to enhanced or reduced 
star formation rates. As shown by various observations, star 
formation in the GC is generally thought to be inefficient, 
i.e. having a factor of > 10 smaller star formation rate per 
free-fall time than what is inferred for local star-forming re¬ 


gions ( Lis et al.|I994[pV[urray|2QII| Kauffmann et al.]|20I3 

Longmore et al.| 2013a). If indeed the high level of turbu¬ 
lence and strong heating by the ISRF and by cosmic rays 
renders GC clouds very inefficient at forming stars, then this 
would provide a simple explanation for the low efficiency of 
star formation in the region as a whole. 

Our simulations demonstrate that the rate and effi¬ 
ciency at which star formation occurs depends on the a 
parameter of the cloud. It is therefore useful to look at the 
typical a values inferred from observations in the GC. It 
turns out that most clouds and clumps in the GC tend to 
be rather unbound, i.e. having values of a > I. For exam¬ 
ple, we find a mean value a = 3.0± 1.6 by computing a for 
different clumps in the GC cloud GCM-0.02-0.07 based on 


data given by Tsuboi & Miyazaki (2012). Kauffmann et al. 


( 20I3| ) evaluated a for the entire GC cloud GO.253+0.016, 
finding a value of a ~ 3.8 + I.O. However, [Longmore et ah] 


(2012) derived a different value. They estimate this cloud 
to be roughly in virial equilibrium with a ^ 1. In this case, 
the difference between their estimate and the IKauffmannI 


et al. (2013) value comes from their decision to exclude an 


additional velocity component in the calculation of a. More 


recently, Rathborne et al. (2015) have computed the virial 


parameter using data from various observed molecular tran¬ 
sitions, finding that the outer regions of GO.253+0.016 may 
be unbound, while its central region may be bound and col¬ 
lapsing (see Fig. 14 in Rathborne et al.|20I5 ). Furthermore, 
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if we evaluate a for numerous cores from data given in Ta¬ 
ble 2 in Johnston et al. (2014) for the Brick, we also find 
that most cores tend to be unbound with a > 2.0. Overall, 
therefore, it seems plausible that many of the clouds in the 
GC that are not currently forming stars have a of a few. 

Regarding measurements of the star formation efficien¬ 
cies per free-fall time, Murray (2011) for example observed 
a Galaxy-wide average value of Cff 0.6%, which is compa¬ 
rable to the lowest in our simulations for highly unbound 
clouds. Generally, a typical estimate for the star formation 
efficiency per free-fall time in the Galaxy is ~ 1%, i.e. star 
formation is quite slow in GMGs ( Murray ||20 11 Krumholz 


et al.|2012 ). We recover larger values than this in almost all 


of our clouds, rather than hnding values strongly suppressed 
compared to the Galactic average, as would be required to 
explain the low efficiency of star formation at the galactic 
center. This holds for runs with both low {a — 0.5) and high 
(a = 8.0) values of a, where Cff > 1% for the two models 
with different initial number density. However, we note that 
these values reflect observed average efficiencies of star for¬ 
mation in the Galaxy as a whole. Thus, a direct comparison 
with the SFEs derived from numerical simulations of clouds 
in isolation presented in this study is complicated, since we 
are only focussing on one specific physical cloud realization, 
instead of a larger number of individual MGs. 


4.2 Comparison to previous studies of unbound 
clouds 

Several previous studies tried to reveal the dependence of 
the SFE on the virial a parameter. Eor example, |Glark et al.| 
(2008) altered the initial level of turbulent support and hnd 


that a wide range of SEEs are possible. The SFEs in their 
study range from up to 60% to as low as around 0.3% after 


two free-fall times. Furthermore, Glark et al. (2008) also 


observed a decreasing SFE with increasing a, as confirmed 
in our study. A similar result was found by [Bonnell et al.| 
(2011), who investigated the formation of young stars in a 


single MG with a total mass that is a factor of 10 smaller 


than our clouds. Bonnell et al. (2011) concluded that even 


small changes in the binding energy of the cloud can cause 
large variations in terms of the SFE, similar to what we find 
in this paper. They gave an overall SFE of 15% when their 


calculation was terminated. Moreover, Glark et al.l (2005) 


also showed that unbound clouds result in inefficient star 


formation, in agreement with 

Klessen et al.| (2000b), 

Heitsch 

et al. (2001), IGlark & Bonnell (2004), Glark et al. 

to 

o 

o 

00 

Bonnell et al.|(2011) and the results in this paper. 


Furthermore, 

Padoan et al. (2012) analyzed the SFR in 


supersonic MHD turbulence, also finding that Cff decreases 
exponentially with increasing a. In particular, they found 
that the SFE is insensitive to changes in the sonic Mach 
number, but sensitive to the Alfvenic Mach number. In their 
study, a decrease in the Alfvenic Mach number (equivalent 
to an increase of the magnetic field) additionally reduced 
the SFE, but only by a factor less than 2. Hence, although 
we do not account for MHD turbulence, we conclude that an 
additional external magnetic field might also slightly reduce 
our values of eff. These findings are also in agreement with 


previous studies by Klessen et al. (2000b) and Heitsch et al. 


(2001), who analyzed the gravitational collapse in turbu¬ 


lent molecular clouds. They find that star formation cannot 


be prevented by MHD turbulence, but the magnetic fields 
delay the local collapse due to the magnetic pressure. Fur¬ 
thermore, they also hnd that strong turbulence can provide 
some support on global scales, but may trigger collapse and 
star formation locally at the stagnation points of convergent 
shocks, in analogy to the hndings in this study. 

All these studies analyzed the impact of the turbulent 
kinetic energy on the SFE. However, regarding the generally 
more extreme physical conditions in the GG, we have to use 
a signihcantly higher initial number density and a stronger 
ISRF/GRF than what was used in these previous studies in 
order to better match the environmental parameters found 
in a typical GMZ-like cloud. Surprisingly, our results are 
similar to those found in previous studies, although we have 
made use of a more extreme physical setup regarding the 
internal velocity, the density and the radiation field. 


4.3 Suppressing star formation 


The cloud GO.253+0.016 is supposed to be a typical MG 


in the GG. However, as shown by Kauffmann et al. (2013) 


and Johnston et al. (2014), GO.253+0.016 has almost no 


evidence of current star formation. In contrast to that, our 
numerical models show active star formation independent of 
the initial number density and the virial a parameter of the 
cloud. Indeed, stars form more rapidly in our model clouds 
than appears to be the case in the Galaxy as a whole. While 
our neglect of stellar feedback probably explains some of this 
discrepancy, stellar feedback can only be effective once star 
formation is ongoing, but cannot explain an almost complete 
lack of star formation in GO.253+0.016. We are therefore 
led to conclude that the harsh GG environment and the 
high level of turbulence present in the GG clouds cannot by 
themselves produce low enough star formation efficiencies 
to explain the globally low efficiency of star formation in 
the GG. They also do not seem to be able to explain why 
GO.253+0.016 and its neighbouring clouds are not currently 
forming stars. What then does suppress star formation in 
the GG? 

One possible explanation could be that the turbulent 
velocity held in many of the GG clouds is composed of 
a different mixture of solenoidal and compressive modes 
than in the model clouds in our study (see, e.g. Federrath 
et al. 2010a, Federrath 2013, Federrath Klessen] 20137^ 
Fundamentally, the reason that high levels of turbulence do 
not completely suppress star formation in our model clouds 
is that the same turbulent motions that support the cloud as 
a whole against collapse also compress some of its gas up to 
high densities. This high density gas is formed at stagnation 
points of the turbulent how and is not supported against 
gravitational collapse. It is therefore able to form stars ef- 
hciently. However, if compressive modes are absent and the 
velocity held is dominated by purely solenoidal turbulence 
(generated e.g. by the strong shear experienced by the clouds 
as they orbit the center of the Galaxy), then less gas will be 
compressed to high densities, and it is plausible that star 
formation could be more strongly suppressed. In addition, 
the strong magnetic held present in the GG might reduce 
the level of star formation by slowing down the collapse of 
dense cores ( Pillai et al.||2015 ). However, as mentioned be¬ 
fore, this typically leads to a reduction in the star formation 
rate and efficiency of only about a factor of ~ 2, see e.g. 
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Peters et al. ( 2011|), Hennebelle et al. (2011), Commergon 


et al. (2011) or Seifried et al. (2013). 


Alternatively, it might be that the idea that we can ex¬ 
plain the low SFE of the GC region as a whole by requiring 
the individual clouds to all have low SFEs is simply incor¬ 
rect. Even in our highly turbulent model clouds, star forma¬ 
tion does not begin at t = 0 - there is a brief period in the 
evolution of the cloud during which no stars are yet forming. 
If, as Longmore et al. (2013b) suggest, GO.253+0.016 and its 
neighbouring clouds have only formed very recently, then it 
may simply be that we are catching them too early in their 
lives to have started forming significant numbers of stars. In 
this case, we would need to look elsewhere for an explana¬ 
tion of the galactic center’s low star formation rate. In this 
context, larger-scale effects such as the orbital dynamics of 
the gas in the deep potential well of the GG, feedback from 
massive stars in the form of winds and supernovae, or the in¬ 
flow/outflow of gas into/out of the GG might play important 
roles by helping to suppress the formation of dense clouds in 
the GG region. We will explore the possible impact of these 
various physical processes in a follow-up study. 


4.4 Limitations of the model 


There are a few limitations inherent to our numerical models 
that one should keep in mind when interpreting our results. 
Most notably, we do not model feedback from the stars (e.g. 
jets, stellar winds or radiation) that form during the individ¬ 
ual runs. This would help to reduce CAt and eff, but cannot 
entirely suppress star formation. Furthermore, we do not 
account for spatial variations in the ISRF or the GRF. 

We also deliberately do not account for other impor¬ 
tant physical effects that might modify the star formation 
rate and efficiency, such as large-scale dynamics (e.g. spiral 
arms or spiral instabilities), magnetic fields, galactic tides 
and shear, supernova feedback or the inflow/outflow of gas. 
This is because we only want to look at the effects of turbu¬ 
lence and a high ISRF/GRF on cloud dynamics in isolation. 
We leave the analysis of simulations with further physical 
effects for future studies. 

Moreover, the high computational cost of Arepo simu¬ 
lations of clouds with mean densities n ^ 10 ^ cm“^ means 
that we have to focus on models with number densities 
no = 100 cm“^ and no = 1000 cm“^. However, the density 
of a typical GG cloud like GO.253+0.016 can be even higher 
than several ~ 10 ^ cm“^, as shown by various studies in the 


past ( 

Lis et al. 2001 |Immer et al. 2012 Longmore et al. 

2012 | 

Kauffmann et al.||2013| |Glark et al.||2013| |Johnston 

et al.|2014). Therefore, as shown in Section | 

a 

, we also run 


finding a SFE of eff = 0.7%. Nevertheless, regarding the 
trends given in Table we expect other high density runs 
with a virial parameter lower than 16.0 to form stars at an 
even higher rate due to a lower amount of turbulent kinetic 
energy. 

Fig. provides support for this assumption. It clearly 
shows that eff depends on the mean density no roughly as 
eff ~ (see Section 3.1). If this trend continues to higher 

density, then typical GG clouds with n ~ 10 ^ cm“^ should 
have about ~ 1/3 lower star formation efficiencies per free- 
fall time than our no = 1000 cm“^ model. In this case, we 
still obtain values of eff > 1 %, which are higher than the 


average SFE in the Milky Way (see Section 4.1). This holds 
at least for all runs except for the high a = 8.0 run, for 
which we would obtain eff ~ 0 . 2 %, smaller than the value 
of eff ~ 0.7% found for the a — 16.0 run. However, it is 
unclear whether this single example can be taken as a sys¬ 
tematic evidence for a breaking of our density trends, given 
the stochasticity with its large uncertainties seen among the 
results for the 100 cm“^ and 1000 cm“^ simulations. 


5 SUMMARY AND CONCLUSIONS 


We have performed numerical simulations of molecular 
clouds with the moving mesh code Arepo ( Springel|| 2010 ) 
using environmental properties comparable to those experi¬ 
enced by typical galactic center (GC) clouds. We adopted 
values for the interstellar radiation field (ISRF) and the 
cosmic ray flux (GRF) that are a factor of ~ 1000 larger 


than the values measured in the solar neighbourhood (Clark 
et al.||2013 ). We simulated clouds with initial number den¬ 
sities of no = 100 cm“^ and 1000 cm“^ using different virial 
a parameters of a = 0.5,1.0, 2.0,4.0 and 8.0 for each den¬ 
sity. The total mass was set to a constant value of Mtot = 
1.3 X 10^ M 0 . In addition, we also ran one more extreme 
simulation with an initial number density of 10 ^ cm“^ and 
a virial parameter of a — 16.0. Furthermore, we also did 
two control runs with no = 1000 cm“^ as well as a = 0.5 
and a = 8.0 and with ISRF and GRF parameters geared 
towards the solar neighbourhood. An overview of our model 
parameters is provided in Table We report the following 
findings: 


• We find active star formation with eff > 1% in all mod¬ 
els regardless of the choice of no and a. 

• Our values are more comparable to Galaxy-wide SFEs 
than to the inferred SFE in the GC, which observations 
suggest is a factor > 10 smaller. Star formation is more 
efficient at lower a values, i.e. when the velocity dispersion 
in the cloud is small. 

• The efficiency of star formation decreases by a factor of 
~ 4 — 10 as we increase the virial parameter from a = 0.5 
to a = 8 . 0 . 

• We fit exponential functions eff oc exp(—ca) to the 

data, finding c ~ 0.20 ± 0.02 and c ~ 0.36 ± 0.07 for the 
no = 100 cm“^ and no = 1000 cm“^ density models, respec¬ 
tively. Thus, we find strong evidence that eff depends not 
only on the virial state of the cloud, but also on its den¬ 
sity. To illustrate this quantitatively, we derive a relation 
eff ~ for the star formation efficiencies as a function 

of the density. 

• For virialized clouds, we find that even a lOOOx higher 
ISRF strength and GRF has only a small effect on the 
star formation efficiency, decreasing it by around 20 % com¬ 
pared to the value we obtain for a similar cloud in a solar 
neighbourhood environment. For highly unbound clouds, the 
stronger ISRF and higher GRF at the galactic center has a 
much greater effect, decreasing the star formation efficiency 
by around a factor of 6 . 

• Even in our most extreme models, we find star forma¬ 
tion efficiencies per free-fall time that are close to 1%. None 
of our models produce values that are consistent with the 
low SFE per free-fall time that is inferred for the galactic 
center region as a whole. 
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We therefore conclude that the idea tested in this paper 
- that the high levels of turbulence present in the GC re¬ 
gion, together with the strong ISRF and high CRF combine 
to yield a persistently low star formation efficiencies within 
the dense clouds in this region - does not appear to work 
in practice. It is possible that including additional physical 
ingredients (such as magnetic fields, stellar feedback or re¬ 
alistic orbital parameters around the GC) could reduce the 
star formation efficiencies within individual dense clouds to 
levels that are consistent with the mean value inferred for 
the GC region as a whole. Alternatively, it could be that 
the idea that we can explain the low SFE of the GC region 
as a whole by requiring the individual clouds to all have 
low SFEs is incorrect, and that the bottleneck for star for¬ 
mation in the region is actually the assembly of the dense 
clouds themselves. This is in agreement with |Kruijssen et al.| 
(2014), who speculate that the rate-limiting factor for star 
formation is the slow evolution of the gas towards collapse. 
Distinguishing between these two possibilities awaits further 
work on this topic. 
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APPENDIX A: SIMULATIONS WITH A 
DIFFERENT RANDOM SEED 


We ran n umer ical simulations with the same setup described 
in Section 2.4 for an initial number density of no = 100 cm“^ 


using a different random seed for the turbulent velocity field. 
Eig. |A1| shows the mass of gas that is converted to stars 
(sink particles) for the different a parameters as a func¬ 
tion of time and their accretion rates. Table |A1| gives an 
overview of the star formation efficiencies and the number 
of sink particles that are formed in each cloud. Since star 




Time [Myr] 

Figure Al. Same as Fig. [^and[^ but using a different random 
seed for the turbulent velocity field. The number of sink particles 
formed by the end of each run, as well as the SFEs, are given in 
Table IaTI 


Model name 

€At 

[%] 

eff 

[%] 

Asink 

f. 

[Myr] 

fend 

[Myr] 

At 

[Myr] 

GC-0.5-100-SEED 

5.7 

24.3 

2483 

1.70 

2.73 

1.03 

GC-1.0-100-SEED 

4.1 

15.0 

1420 

1.98 

3.10 

1.12 

GC-2.0-100-SEED 

9.4 

20.7 

2717 

2.40 

4.40 

2.00 

GG-4.0-100-SEED 

7.6 

10.0 

2067 

3.20 

6.56 

3.36 

GG-8.0-100-SEED 

1.1 

1.8 

405 

4.25 

7.00 

2.75 


Table Al. Same as Table but using a different random seed 
for the turbulent velocity field. 


formation sets in later in this turbulent environment com¬ 
pared to our fiducial models presented in Section |3.1| we 
let the two runs GC-4.0-100-SEED and GC-8.0-100-SEED 
evolve beyond one free-fall time in order to get a significant 
number of sink particles. Although star formation starts at 
later times, caused by the different statistical flows of the 
turbulent velocity field, we again measure large efficiencies 
Cff for all numerical models even in highly turbulent environ¬ 
ments. We also find the same statistical trends as already 
observed in Table i.e. that the SEEs per free-fall time 
strongly depend on the virial parameter. However, a direct 
comparison of the Cff values in models with various random 
seeds is complicated due to the statistics of the turbulent ve¬ 
locity fields, leading to a completely different star formation 
history and thus to variable SEE. 


© 0000 RAS, MNRAS 000 , 000-000 



















15 


REFERENCES 

Ballesteros-Paredes J., Klessen R. S., Mac Low M.-M., 
Vazquez-Semadeni E., 2007, Protostars and Planets V, 
pp 63-80 

Bate M. R., Bonnell L A., Price N. M., 1995, MNRAS, 277, 
362 

Bonnell L A., Smith R. J., Clark P. C., Bate M. R., 2011, 
MNRAS, 410, 2339 

Clark R C., Bonnell L A., 2004, MNRAS, 347, L36 
Clark R C., Bonnell L A., Klessen R. S., 2008, MNRAS, 
386, 3 

Clark P. C., Bonnell L A., Zinnecker H., Bate M. R., 2005, 
MNRAS, 359, 809 

Clark R C., Glover S. C. O, Klessen R. S., 2012, MNRAS, 
420, 745 

Clark P. C., Glover S. C. O., Ragan S. E., Shetty R., 
Klessen R. S., 2013, ApJ, 768, L34 
Commercon B., Hennebelle P., Henning T., 2011, ApJ, 742, 
L9 

Draine B. T., 1978, ApJS, 36, 595 
Elmegreen B. G., Scalo J., 2004, ARAVA, 42, 211 
Eederrath C., 2013, MNRAS, 436, 1245 
Eederrath C., Klessen R. S., 2013, ApJ, 763, 51 
Eederrath C., Roman-Duval J., Klessen R. S., Schmidt W., 
Mac Low M.-M., 2010a, AVA, 512, A81 
Glover S. C. O., Clark P. C., 2012a, MNRAS, 421, 9 
Glover S. C. O., Clark P. C., 2012b, MNRAS, 421, 116 
Glover S. C. O., Eederrath C., Mac Low M.-M., Klessen 
R. S., 2010, MNRAS, 404, 2 

Glover S. C. O., Mac Low M.-M., 2007, ApJS, 169, 239 
Greif T. H., Springel V., White S. D. M., Glover S. C. O., 
Clark P. C., Smith R. J., Klessen R. S., Bromm V., 2011, 
ApJ, 737, 75 

Glisten R., Walmsley C. M., Pauls T., 1981, AVA, 103, 197 
Habing H. J., 1968, Bull. Astron. Inst. Netherlands, 19, 421 
Heitsch E., Mac Low M.-M., Klessen R. S., 2001, ApJ, 547, 
280 

Hennebelle P., Commercon B., Joos M., Klessen R. S., 
Krumholz M., Tan J. C., Teyssier R., 2011, AVA, 528, 
A72 

Immer K., Menten K. M., Schuller E., Lis D. C., 2012, 
AVA, 548, A120 

Jappsen A.-K., Klessen R. S., Larson R. B., Li Y., Mac 
Low M.-M., 2005, AVA, 435, 611 
Johnston K. G., Beuther H., Linz H., Schmiedeke A., Ra¬ 
gan S. E., Henning T., 2014, AVA, 568, A56 
Kauffmann J., Pillai T., Zhang Q., 2013, ApJ, 765, L35 
Kennicutt R. C., Evans N. J., 2012, ARAVA, 50, 531 
Kennicutt Jr. R. C., 1998, ApJ, 498, 541 
Klessen R. S., 2000a, ApJ, 535, 869 

Klessen R. S., Glover S. C. O., 2014, ArXiv e-prints, 
1412.5182 

Klessen R. S., Heitsch E., Mac Low M.-M., 2000b, ApJ, 
535, 887 

Kritsuk A. G., Norman M. L., Wagner R., 2011, ApJ, 727, 
L20 

Kruijssen J. M. D., Dale J. E., Longmore S. N., 2015, MN¬ 
RAS, 447, 1059 

Kruijssen J. M. D., Longmore S. N., Elmegreen B. G., Mur¬ 
ray N., Bally J., Testi L., Kennicutt R. C., 2014, MNRAS, 
440, 3370 


Krumholz M. R., Dekel A., McKee C. E., 2012, ApJ, 745, 
69 

Krumholz M. R., Kruijssen J. M. D., 2015, ArXiv e-prints, 
1505.07111 

Krumholz M. R., McKee C. E., 2005, ApJ, 630, 250 
Krumholz M. R., Tan J. C., 2007, ApJ, 654, 304 
Lis D. C., Menten K. M., 1998, ApJ, 507, 794 
Lis D. C., Menten K. M., Serabyn E., Zylka R., 1994, ApJ, 
423, L39 

Lis D. C., Serabyn E., Zylka R., Li Y., 2001, ApJ, 550, 761 
Longmore S. N., Bally J., Testi L., Purcell C. R., Walsh 
A. J., Bressert E., Pestalozzi M., Molinari S., Ott J., 
Cortese L., Battersby C., Murray N., Lee E., Kruijssen 
J. M. D., Schisano E., Elia D., 2013a, MNRAS, 429, 987 
Longmore S. N., Kruijssen J. M. D., Bally J., Ott J., Testi 
L., Rathborne J., Bastian N., Bressert E., Molinari S., 
Battersby C., Walsh A. J., 2013b, MNRAS, 433, L15 
Longmore S. N., Rathborne J., Bastian N., Alves J., As- 
censo J., Bally J., Testi L., Longmore A., Battersby C., 
Bressert E., Purcell C., Walsh A., Jackson J., Poster J., 
Molinari S., Meingast S., Amorim A., Lima J., 2012, ApJ, 
746, 117 

Mac Low M.-M., Klessen R. S., 2004, Reviews of Modern 
Physics, 76, 125 

McKee C. E., Ostriker E. C., 2007, ARAVA, 45, 565 
Molinari S., Bally J., Noriega-Crespo A., Compiegne M., 
Bernard J. P., Paradis D., Martin P., Testi L., Barlow M., 
Moore T., Plume R., Swinyard B., Zavagno A., Calzoletti 
L., Di Giorgio A. M., Elia 2011, ApJ, 735, L33 
Murray N., 2011, ApJ, 729, 133 
Nelson R. P., Langer W. D., 1997, ApJ, 482, 796 
Padoan P., Haugbplle T., Nordlund A., 2012, ApJ, 759, L27 
Peters T., Banerjee R., Klessen R. S., Mac Low M.-M., 
2011, ApJ, 729, 72 

Pillai T., Kauffmann J., Tan J. C., Goldsmith P. E., Carey 
S. J., Menten K. M., 2015, ApJ, 799, 74 
Rathborne J. M., Longmore S. N., Jackson J. M., Poster 
J. B., Contreras Y., Garay G., Testi L., Alves J. E., Bally 
J., Bastian N., Kruijssen J. M. D., Bressert E., 2014, ApJ, 
786, 140 

Rathborne J. M., Longmore S. N., Jackson J. M., Kruijssen 
J. M. D., Alves J. E., Bally J., Bastian N., Contreras Y., 
Poster J. B., Garay G., Testi L., Walsh A. J., 2015, ArXiv 
e-prints, arXiv: 1501.07368 

Scalo J., Elmegreen B. G., 2004, ARAVA, 42, 275 
Schmidt M., 1959, Astrophys.J., 129, 243 
Schneider N., Andre P., Konyves V., Bontemps S., Motte 
E., Eederrath C., Ward-Thompson D., Arzoumanian D., 
Benedettini M., Bressert E., Didelon P., Di Erancesco J., 
Griffin M., Hennemann M., Hill T., 2013, ApJ, 766, L17 
Seifried D., Banerjee R., Pudritz R. E., Klessen R. S., 2013, 
MNRAS, 432, 3320 

Sembach K. R., Howk J. C., Ryans R. S. L, Keenan E. P., 
2000, ApJ, 528, 310 

Shetty R., Kelly B. C., Bigiel E., 2013, MNRAS, 430, 288 
Shetty R., Kelly B. C., Rahman N., Bigiel E., Bolatto A. D., 
Clark P. C., Klessen R. S., Konstandin L. K., 2014, MN¬ 
RAS, 437, L61 

Smith R. J., Glover S. C. O., Clark P. C., Klessen R. S., 
Springel V., 2014a, MNRAS, 441, 1628 
Smith R. J., Glover S. C. O., Klessen R. S., 2014b, MNRAS, 
445, 2900 


© 0000 RAS, MNRAS 000 , 000-000 



16 Bertram et al. 


Springel V., 2010, MNRAS, 401, 791 

Tmelove J. K., Klein R. L, McKee C. F., Holliman II J. H., 
Howell L. H., Greenough J. A., Woods D. T., 1998, ApJ, 
495, 821 

Tsuboi M., Miyazaki A., 2012, Publications of the Astro¬ 
nomical Society of Japan, 64, 111 
Yusef-Zadeh F., Wardle M., Roy S., 2007, ApJ, 665, L123 


© 0000 RAS, MNRAS 000 , 000-000 



